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We report the results of photometric monitoring of the Melotte 66 field in BVI 
filters. Ten variables were identified with nine being new discoveries. The sample in- 
cludes eight eclipsing binaries of which four are W UMa type stars, one star is a 
candidate blue straggler. All four contact binaries are likely members of the cluster 
based on their estimated distances. Ten blue stars with U — B < —0.3 were detected 
inside a 14.8 x 22.8 arcmin^ field centred on the cluster. Time series photometry for 
7 of them showed no evidence for any variability. The brightest object in the sample 
of blue stars is a promising candidate for a hot subdwarf belonging to the cluster. 
We show that the anomalously wide main sequence of the cluster, reported in some 
earlier studies, results from a combination of two effects: variable reddening occuring 
across the cluster field and the presence of a rich population of binary stars in the 
cluster itself. The density profile of the cluster field is derived and the total number 
of member stars with 16 < 1/ < 21 or 2.8 < My < 7.8 is estimated conservatively at 
about 1100. 
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1 INTRODUCTION 

Melotte 66 (a = 07''26™21.'9, 5 = -47°41'19", J2000) be- 
longs to a small sample of old galactic open clusters (Kassis 
et al. 1997). It is located at a relatively large galactic lat- 
itude (/ = 259?6, b = — 14?3). A pioneering study, based 
on the photographic photometry was conducted by Eggen 
& Stoy (1962). Subsequent investigations based on pho- 
tographic and photoelectric photometry (Hawarden 1976; 
Anthony- Twarog et al. 1979) were followed by several pa- 
pers published in the last two decades of the 20th century. 
Kassis et al. (1997) used deep VI photometry to derive an 
age of 4±1 Gyr and distance modulus {m — M)o ~ 13.2^0 '^. 
Anthony- Twarog et al. (1994) used vbyHjS data to discuss 
the possible causes of an atypically wide cluster main se- 
quence which was originally noted by Anthony- Twarog et 
al. (1979). They excluded differential reddening across the 
cluster field as a cause of this effect. 

The cluster stands out from the sample of the known 
old open clusters in its exceptionally low metallicity. Twarog 
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et al. (1995) derived [Fe/H] = -0.53 ± 0.08 based on UBV 
photometry of turnoff stars. Friel & Janes (1993) obtained 
[Fe/H] = —0.51 ± 0.11 from a spectroscopic analysis of 4 
giants while Gratton & Contarini (1994) derived [Fe/H] = 
—0.38 ± 0.15 from high resolution spectra of two giants. 

This paper is a contribution to the systematic search 
for short period variables in open clusters conducted by our 
group over the last two decades. One of the goals is to es- 
tablish a relation between age and a relative frequency of 
occurrence of contact binaries in stellar clusters. A sum- 
mary of some of our earlier results can be found in Rucinski 
(1998). Here we report the results of a survey for variable 
stars in the field of Melotte 66 and also present deep CCD 
UBV I photometry for the cluster. 



2 OBSERVATIONS AND REDUCTIONS 

The field of Melotte 66 was surveyed for variable stars with 
the 1-m Swope telescope at Las Campanas Observatory. The 
observations were collected on a total of 21 nights during 
4 observing runs during the period February-March 1992. 
Two different cameras were used: a 1024 x 1024 Tektronics 
CCD with a scale of 0.61 arcsec pixel"^ (TEK2 camera with 
a 10.4 X 10.4 arcmin^ field) and a 2048x2048 Ford Aerospace 
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Table 1 . Summary of Melotte 66 observations 



Run 


Dates 


Nights 


No of 


exposures 


Median seeing 






no 


B 


V 


I 


in V filter 


FORD2-1 


03-08.02.1992 


5 


8 


59 




1"5 


FORD2-2 


20-24.03.1992 


3 




12 


44 


1"2 


TEK2-1 


09-16.02.1992 


7 




134 


5 


1"5 


TEK2-2 


13-18.03.1992 


6 


68 


66 




1"5 



CCD with a scale 0.435 arcsec pixel" ^ (FORD2 camera with 
a 14.8 X 14.8 arcmin^ field). 

Most of the images were collected with a V filter with 
exposure times ranging from 60 s to 480 s with a median 
value of 420 s. Exposures were also obtained in the B and 
I bands. A summary log of the observations is listed in Ta- 
ble 1. An additional set of UBVI observations to calibrate 
the photometry and to construct a color-magnitude diagram 
(CMD) for the cluster field was obtained on the nights of 
1999 November 17-20 (UT) using the 2048 x 3150 SITE3 
camera. With a scale of 0.435 arcsec pixel" ^ this camera pro- 
vides a field of view of 14.8 x 22.8 arcmin^ . Several exposures 
of different length were obtained in each of the four filters. 
All images were corrected for the known non-linearity of the 
SITE3 camera using the procedure described in Hamuy et 
al. (2006). 

Preliminary processing of the CCD frames was 
done with standard routines in the IRAF-CCDPROCQ 
package. Profile photometry was extracted using the 
DAOPHOT/ALLSTAR package (Stetson 1987). For each 
camera/filter combination a "master" frame was selected. 
This was used to create a reference list of objects to be mea- 
sured in the fixed-position mode on the remaining frames. 
Instrumental magnitudes were transformed to the system 
defined by the master frame. The resulting data bases were 
searched for variable objects using codes using the analy- 
sis of variance (AoV) statistic (Schwarzenberg-Czerny 1996) 
and the AOVTRANS algorithm (Schwarzenberg-Czerny & 
Beaulieu 2006). AoV periodograms were calculated for pe- 
riods spanning the range from 0.05 days to 20 days. After 
the rejection of some spurious detections we ended up with 
a list of 10 certain variable stars. Nine of these are new dis- 
coveries while one is an eclipsing binary V345 Pup reported 
by Kaluzny and Shara (1988). Equatorial coordinates of the 
variables were determined using 1651 stars from the USNO 
A-2 catalog (Monet et al. 1998) which were identified on the 
l/-band image obtained with the SITE3 camera. These co- 
ordinates are listed in Table 2 along with the coordinates of 
a UV bright star labelled CI (see Sec. 3.1) The last column 
of Table 2 gives the angular distance from the cluster cen- 
tre as determined from our data (RA(2000) = 07''26"'21.=9 
and DEC(2000) = -47°41'19"). In Fig.[T]we present finding 
charts for the stars listed in Table 2. 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under a cooperative agreement with 
the National Science Foundation. 



Table 2. Equatorial coordinates of variables and a UV bright 
star in the field of Melotte 66 



ID 


"2000 [deg] 


52000 [deg] 


r [arcmin] 


VI 


111.59873 


-47.70060 


0.78 


V2 


111.53339 


-47.62291 


4.58 


V3 


111.61123 


-47.67409 


1.18 


V4 


111.57849 


-47.69664 


0.71 


V5 


111.56815 


-47.64329 


2.87 


V6 


111.72366 


-47.76692 


7.12 


V7 


111.56230 


-47.72299 


2.38 


V8 


111.54562 


-47.71756 


2.54 


V9 


111.71073 


-47.75736 


6.35 


VIO 


111.63718 


-47.68840 


1.85 


CI 


111.57537 


-47.74433 


3.31 



2.1 Photometric calibration 

We have used the data collected on the night of 1999 Novem- 
ber 18 to calibrate our photometry. Observations of the clus- 
ter field were bracketed by observations of three Landolt 
fields containing a total of 28 standard stars with BVI mag- 
nitudes (Landolt 1992; Stetson 2000) and 19 stars with U 
magnitudes (Landolt 1992). The T Phe field was observed 
twice while each of the RU 149 and RU 152 fields was ob- 
served once. Standards were observed at air-masses spanning 
the range 1.12-1.24 while Melotte 66 was observed at an air- 
mass of about 1.1. Average extinction coefficients for Las 
Campanas were assumed in determining the linear transfor- 
mations between the instrumental and the standard system. 
The total uncertainties of the zero points of our photometry 
are about 0.02 mag for BVI filters and about 0.1 mag for 
the U filter. These uncertainties include errors in the aper- 
ture corrections derived for the frames of Melotte 66. The 
relatively poor quality of the [/-band transformation can be 
explained by large differences in the UV spectral responses 
of the SITe3 CCD camera and the RCA 1P21 photomulti- 
plier defining the UBV system (Johnson 1963). The median 
values of the formal internal errors of our photometry for 
the cluster range from av = 0.010 and av-i = 0.021 at 
V = 15.5 to av = 0.018 and av-i = 0.027 ai V = 19.5. 

We compare our VI photometry with the data taken 
from Kassis et al. (1997) in Fig. [5] The mean difference (in 
the sense of "our" measurements minus "theirs") is — 0.005± 
0.001 and -0.023 ± 0.001 for V a.nd V - I, respectively. 
There is no clear color dependence of the residuals for the V 
magnitudes. There are some systematic trends for the V — I 
residuals. 
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Figure 1. Finding charts for variables VI— VIO and for the UV bright star CI. Each chart is 1 arcmin on a side, with east to the left 
and north up. 
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Figure 2. Residuals of Y and V — 1 for this work and Kassis et 
al. (1997) as a function of V — /. 



3 RESULTS FOR VARIABLES 

Table 3 lists some basic characteristics of the light curves of 
variables Vl-VlO. For each star we list V magnitude and col- 
ors measured at maximum light. The full range of observed 
magnitudes in the V band is listed as AV^. The location of 
the variables on the cluster V/B — V and V jV — I diagrams 
is shown in Fig. |3] while their light curves are presented in 
Figs. I3HS1 Stars V3-10 show periodic variability and can be 
classified as eclipsing binaries (EC). Moreover, stars V4, V6, 
V8 and V9 are likely contact binaries (EW) based on their 



V 
14 



15 
16 
17 
18 
19 
20 



VIA 



AV2 



;-*V4;., 



■ye. 



V8a - 



VIO— ^ > 




0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 1.2 
V-l B-V 

Figure 3. Color-magnitude diagrams of Melotte 66 field with 
marked positions of variables. 



periods, colors and light curves. Variables VI and V2 are 
located on or near the cluster subgiant branch while eclips- 
ing binary V3 is a candidate blue straggler. Three other 
non-contact binaries are located on or slightly above the 
cluster main sequence. Unfortunately, none of these seem 
to be a good candidate for spectroscopic follow-up aimed 
at a determination of cluster distance and age. V5 shows 
shallow eclipses while VIO is too faint for high resolution 
spectroscopy. The light curve of V7 is well defined and its 
shape indicates that the system is a likely semi-detached 
algol. 

All variables are located inside the cluster radius which 
we estimate at 9. '2 in Sec. 4. An examination of the angular 
distances listed in Table 2 shows that most of the variables 
are in the central part of Melotte 66. There is radial velocity 
information available for only one of the variables. Collier & 
Reid (1987) list Vr = 75 ± 7 km s"^ for V2=#1210 based on 
a single observation. The mean value for 59 stars considered 
to be cluster members is 44 ± 12 km s~^. This seems to 
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Figure 4. BVI phased light curves of periodic variables V3, V4, V5, V7, V8 and VIO. 
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Figure 5. Phased light curves of periodic variables V6 and V9. 
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Figure 6. Unphascd light curves of variables VI and V2. 



indicate that V2 has a low probability to be a radial velocity 
member of Melotte 66. However, this evidence is weak as the 
radial velocity of V2 is likely to be variable. 

Some indication of membership status can be provided 
for 4 contact binaries. We have applied the absolute bright- 
ness calibration established by Rucinski (2000) to estimate 
their absolute magnitudes. The calibration gives Mv as a 
function of unreddened color and orbital period. We adopted 
reddening of = 0.22 resulting from E{B-V) = 0.16 

advocated by Anthony- Twarog, Twarog & Sheeran (1994). 
Using dereddened V — I colors and periods listed in Table 
3 we obtained Afy = 3.51, My = 3.07, Mv = 3.79 and 
Mv = 6.33 for V4, V6, V8 and V9, respectively. The appar- 
ent distance moduli of these four variables follow from the 
observed values of V-max- We obtained (M — m)v = 13.14, 
(M-m)v = 14.17, {M~m)v = 14.36 and (M-m)v = 13.78 
for V4, V6, V8 and V9, respectively. These values can be 
compared with the apparent distance modulus of the clus- 
ter of (M — m)v = 13.75 as measured by Kassis et al. (1997). 
The formal error of My obtained from Rucinski's calibration 
is about 0.3 mag and so the estimated distance moduli are 
consistent with cluster membership for all four contact bi- 
naries. This conclusion holds if we adopt E{B — V) — 0.23 
for the cluster reddening as implied by maps of Schlegel et 
al. (1998). 

3.1 Search for cataclysmic variables 

So far there are only 3 confirmed and one candidate cata- 
clysmic variables known in the whole sample of galactic open 
clusters: one in M 67 (Gilhland et al. 1991), two in NGC 6791 
(Kaluzny et al. 1997) and one in NGC 2158 (Mochejska et 
al. 2006). All of these clusters are old, and as Melotte 66 is 
an old and rich cluster we decided to search it for possible 
cataclysmic variables. Two methods were used. 



First we used the ISIS image subtraction package 
(Alard & Lupton 1998; Alard 2000) to look for objects show- 
ing outbursts. Cataclysmic variables of dwarf novae type 
have average A/y = 7.5 at quiescence with outburst mag- 
nitudes spanning the range 2-8 mag (Warner 1995). At the 
cluster distance they would be observed at V ~ 21 at min- 
imum light and at 13 < 1^ < 19 at maximum light. With a 
limiting magnitude of our observations of V ~ 21 it should 
be possible to detect objects of this type while they are in 
outburst. Our search gave a negative result. 



The second method used relies on the fact that cata- 
clysmic variables have blue U — B colors. In Fig. [7] we show 
a. V/U ~ B diagram for the cluster field. It contains 10 ob- 
jects with U — B < —0.3. For 7 of these objects we have 
time series photometry obtained with the FORD2 camera 
and for 2 photometry obtained with TEK2 camera. Exami- 
nation of these light curves showed than none of the objects 
shows any convincing evidence for variability. The brightest 
of the blue objects, which we denote CI, has V = 18.57, 
B-V = 0.07 and (7 - B = -0.87. The lack of evidence for 
variability, including comparison of photometry from 1992 
and 1999 seasons, indicates that it is unlikely to be a quasar. 
It is possible that the star is a hot subdwarf belonging to 
the cluster. We note that CI is located very close to the 
projected cluster centre. Very few confirmed hot subdwarfs 
are known in open clusters and it is of interest to obtain a 
spectrum of CI to clarify its nature. If CI is a member of 
the cluster then its absolute magnitude is My = 5.4. This 
value corresponds to the faint end of the absolute magnitude 
distribution observed for hot subdwarfs in the field and in 
globular clusters (Lisker et al. 2005; Moehler et al. 2002, 
2004). 
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Table 3. Parameters of variables from the field of Melotte 66 



ID 






(B — V)iaax 


Ay 


P [days] 


To HJD 2448000+ 


Remarks 


VI 


14.626 


1.29 


1.10 


0.15: 


P ~ 8^* 




periodic? 


V2 


15.884 


1.18 


0.99 


0.1: 








V3 


16.341 


0.45 


0.37 


0.30 


0.8015(2) 


657.7128 


Eel 


V4 


16.65 


0.79 


0.61 


0.17 


0.4020(17) 


656.6764 


Ecl-EW 


V5 


17.03 


0.85 


0.65 


0.15 


0.7413(2) 


658.7648 


Eel 


V6 


17.24 


0.96 




0.12 


0.6974(3) 


657.6511 


Ecl-EW 


V7 


17.519 


0.86 


0.64 


1.78 


0.5942(1) 


656.5508 


Ecl=V345 Pup 


V8 


18.16 


0.71 


0.55 


0.16 


0.32903(6) 


655.5672 


Ecl-EW 


V9 


20.11 


1.29 




0.49 


0.2386(4) 


656.7105 


Ecl-EW 


VIO 


20.2 


1.1 


1.0 


0.9 


0.8882(9) 


668.7426 


Eel 
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Figure 8. Star density as a function of radial distance from the 
center of Melotte 66. 
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Figure 7. U — B versus V color-magnitude diagram of Melotte 
66. 



4 ANALYSIS OF THE CMD 

Some early photometric studies of Melotte 66 revealed an 
unexpectedly large widths for the cluster subgiant branch 
and upper main sequence (Hawarden 1976; Anthony-Twarog 
et al. 1979). Anthony-Twarog et al. (1994) used vbyH/3 CCD 
photometry to eliminate differential reddening and varia- 
tions in cluster metallicity as possible causes of these large 
widths. They suggested a broad range of rotational velocities 
among cluster stars as an explanation. 

To further investigate this issue we started with a de- 
termination of the stellar density profile for the field ob- 
served with the SITES camera. We observed 2134 stars with 
16 < V < 21 located near the cluster main sequence in the 
V/V — I CMD. The resulting density profile is shown in 
Fig. [S] The projected density distribution fiattens at a ra- 
dius of about 550 arcsec. The surface density of the field stars 
can be estimated at 8.24 ± O.lli? — 4 stars per arcsec^. This 
implies that the SITE3 camera observations fail to cover 
the whole cluster E-W on the sky. By integrating the den- 
sity profile shown in Fig. |8] and subtracting the contribution 
from the field stars we conclude that the cluster contains 
1122 stars with 16 < V < 21 or 2.8 < Afv < 7.8 inside a 
radius of i? = 550 arcsec. Note that this is a conservative 



lower limit because no corrections for the incompleteness of 
the photometry have been applied. At a radius of 225 arcsec 
from the cluster centre the cluster surface density is still a 
factor of 3 higher than the field star density. In Fig. |9] we 
show V/V — I CMDs for two groups of stars: those lying 
inside a radius of 7? = 225 arcsec from the cluster centre 
and those lying in the outer part of the observed field at 
R > 550 arcsec. The inner circle and the outer region cover 
equal areas on the sky. For each star from the outer field a 
nearest match in the CMD for the inner ring was located. 
Subsequently a pair of stars with the lowest separation was 
removed from both corresponding lists. This procedure was 
continued until it was impossible to locate pairs with a sep- 
aration SV < 0.25 and S(V - I) < 0.15. 

The resulting "cleaned" CMD for the inner region of 
Melotte 66 is shown in Fig. 1101 One may notice 7 candidate 
blue stragglers as well as a clump of yellow stragglers at 
V ~ 15.4 and V" — / « 0.8. The occurence of such objects 
in the cluster was first noted by Hawarden (1976) and sub- 
sequently discussed by other investigators whose papers are 
quoted above. Fig. [10] shows that the cluster main sequence 
is very sharply defined on the blue side. In particular, for 
17.4 < V < 19.5 it is possible to distinguish a narrow, well 
defined main sequence corresponding to a sample dominated 
by single stars. Above this sequence is a second sequence 
consistent with a sample of binary stars with mass ratios 
close to unity. It merges with the "single" sequence near the 
turn-off region. As can be seen in Fig. [31 the eclipsing bina- 
ries V4, V5, V7 and V9 are located on the binary sequence. 

The presence of a well populated binary sequence in the 
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Figure 9. The CMD's for the field covering central part of 
Melotte 66 (left) and for the "outer" field (right). 
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Figure 10. Field-star corrected CMD for the central part of 
Melotte 66. 



CMD of Melotte 66 has already been noticed by Kassis et 
al. (1997). 

We have estimated the widths of these " single" and " bi- 
nary" sequences using stars from the cleaned CMD in a rect- 
angle given by 17.5 <V < 19.5 and 0.7 < 1^ - / < 1.1. The 
blue edge of the main sequence in this rectangle was fitted by 
second order polynomial. For each star we have calculated its 
distance (in color) from this polynomial. In Fig.[TT]we show 
a histogram of star counts in bins of 0.015 in V — I. Fig. 1111 
shows two peaks representing the two sequences. A double- 
Gaussian fit to the two peaks indicates widths of V — I = 
0.042±0.002 and 0.048±0.005 for the "single" and "binary" 
sequences, respectively. The two sequences are well resolved 
in over the range 17.5 < V < 19.5. At fainter magnitudes 
the sequences are smeared by the uncertainties of the color 
measurements. At 16.0 < V < 17.5 the sequences cross each 
other leading to an apparent broadening of the upper main 
sequence. Fig. 1101 shows that the cluster subgiant branch is 
well populated, with 4 or 5 stars scattered on the blue side. 
We speculate that these stars represent an extension of the 
binary sequence in the turn-off region. 

We have checked how the E{B — V) extinction varies 
across the field of the cluster based on the Schlegel et al. 
(1998) reddening maps. In Fig. 1121 we show an E{B — V) 
differential reddening map with 5 arcmin resolution in a 25 
X 25 arcmin box centred on the cluster. It is clear that even 
inside the cluster radius AE{B — V) can be as large as 0.09. 
According to Kassis et al. (1997) Melotte 66 is located about 
1.1 kpc below the Galactic disk and 4.4 kpc from the Sun. 
The line of sight to the cluster does not cross any outer spi- 
ral arm of the Milky Way and so the observed gradients of 
E(B — V) occur in the interstellar matter located between 
the Sun and the clusterQ Our data indicate the large width 
of the cluster main-sequence can be explained by a high rel- 



A possible extension of the Perseus spiral arm (Caswell & 
Haynes 1987) is present at the Galactic longitude of the cluster 
at a distance of about 5 kpc. However, it is unlikely that this 
contributes to the reddening in the direction of Melotte 66 due 
to the large z of the cluster. 



FWHMgg = 0.042 +/- 0.002 
FWHMeg = 0.048 +/- 0.005 




0.25 



Figure 11. Distribution of A{V — I) distances from the blue 
edge of the cluster main sequence (details in the text). SS and 
BS denote a "single" and "binary" stars sequence respectively. 
Solid curve is a fit to a double-gaussian. Broken curves correspond 
to distinct two gaussians plotted with parameters derived in the 
double-gaussian fit. 



ative frequency of binaries among the cluster stars together 
with differential reddening across the cluster field. 



5 CONCLUSIONS 

We have detected a total of 10 photometric variables in the 
field of the old open cluster Melotte 66. Four out of eight 
eclipsing binaries are contact systems which are probable 
members of the cluster. Like other old open clusters Melotte 
66 seems to have a rather high relative frequency of contact 
binaries. Our estimate of the total number of cluster mem- 
bers with 16 < < 21 at 1122 stars leads to a relative 
frequency of EW stars of 0.36±0.18%. Note, however, that 
the sensitivity of our survey for variables is rather low for 
V > 20. This relative frequency can be compared to that ob- 
served for field stars in the solar vicinity which is estimated 
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Figure 12. Differential E{B — V) reddening map with 5 arcmin 
resolution in 25 x 25 arcmin box centred at the Melotte 66 coor- 
dinates. Circle denotes cluster size R = 550 arcsec. 

at 0.2% for 3.5 < Mv < 5.5 and at 0.1% for My = 6.0 
(Rucinski 2006). Three other binaries, including a candidate 
blue straggler, are detached systems and one is probably a 
semi-detached binary. We also identified a promising candi- 
date for a hot subdwarf cluster member. 

The projected stellar density profile was obtained for 
the cluster field. We show that the angular radius of Melotte 
66 reaches at least 550 arcsec. When corrected for the con- 
tamination by field stars the cluster CMD shows a well de- 
fined and narrow main sequence accompanied by a rich se- 
quence of binary stars. 
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